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Abstract 

We investigate a lattice scalar field theory in the presence of a bias favouring 
the establishment of an energy current, as a model for stationary nonequilibrium 
processes at low temperature in a non-integrable system. There is a transition at 
a finite value of the bias to a gapless modulated phase which carries a classical 
current; however, unlike in similar, integrable, models, quantum effects also allow 
for a non-zero current at arbitrarily small bias. The transition is second order in 
the magnetically disordered phase, but is pre-empted by a first-order transition in 
the ferromagnetic case, at least at the mean-field level. 

1 Introduction 

Finding a general framework for the description of nonequilibrium processes has been a 
long standing problem in statistical physics. In particular, there is a substantial body of 
literature dealing with the microscopic theory of energy, or heat, conduction |J fU @ H ■ 
Unfortunately, no clear conclusions have been reached as a result of these investigations. 
Results are dependent on boundary conditions and on basic assumptions concerning the 
underlying process responsible for energy conduction. 

Recently, Antal, Racz, and Sasvari have proposed studying a related but slightly 
different problem, namely the microscopic description of energy conduction in quantum 
systems at zero temperature |J || . Denoting the quantum hamiltonian by H ^ and the 
space integral of the energy current by J, they introduce a Lagrange multiplier A and 
study the ground state properties of — XJ, for certain exactly solvable one dimen- 
sional quantum spin systems. An interesting feature of these models is that the energy 
current turns on only at a critical, nonvanishing, value of the Lagrange multiplier. At the 
critical point the system undergoes a second order phase transition to a gapless energy- 
conducting state, in which the spin correlation function has a characteristic oscillatory 
behaviour, whose amplitude is modulated by a power law fall-off. 

It is of considerable interest to find out to what extent the phenomena found in []5J , H 
apply to other systems, in particular whether a phase transition to an energy-conducting 
state persists in higher dimensions, and whether integrability plays an essential role in 
the results. 
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Having departed from exactly solvable theories, one needs a systematic approximate 
method to attack the problem. In quantum field theory the effective action 0], com- 
bined with the loop expansion, has been extremely successfully in treating all kinds of 
problems related to symmetry breaking. Standard applications of the method of effec- 
tive action are usually restricted to the breaking of global or gauge symmetries ||. In 
spontaneous symmetry breaking of continuous symmetries in three or more dimensions 
a scalar field acquires a constant vacuum expectation value. As this expectation value, 
or classical field, is constant the effective action reduces to the effective potential, simpli- 
fying calculations considerably. An additional bonus of working with the effective action 
is that the inverse correlation function is obtained as the second functional derivative of 
the effective action with respect to the classical field. The definiteness of the spectral 
function of the propagator implies the convexity of the effective action. Just like its 
analogue in statistical mechanics, the Gibbs free energy, the absolute minimum of the 
effective action selects the correct phase of the system. Although, as we shall show, the 
method of the effective action when applied to our problem leads to a coordinate depen- 
dent classical field and, thus, spontaneous breaking of translation invariance, we expect 
that the effective action is still minimised by the correct classical field configuration. 

The model we consider in this paper is a continuous spin version, in D space di- 
mensions, of the transverse Ising model discussed in for D = 1. In the absence 
of a current bias, the model has a phase transition from a disordered state to one in 
which in which the Z2 symmetry is spontaneously broken, in the universality class of 
the (D + l)-dimensional classical Ising model. The phase diagram in the presence of the 
current bias A is shown schematically in Fig. [I]. Qualitatively it is quite similar to that 
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Figure 1: Schematic phase diagram at D > 1. T signifies translation symmetry. The 
dashed line represents a first order transition. 

found in ||. However, there are important differences. At a finite value of A there is a 
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transition to a situation in which, for D > 1, the ground state is no longer invariant un- 
der lattice translations. In this phase there is a characteristic wave number p Q such that 
Fourier components of the field which are an integer multiple of this acquire a vacuum 
expectation value. For /i 2 > 0, only the fundamental frequency is important near the 
transition, and its amplitude vanishes like (A — Ac) 1 / 2 . This classical spin density wave 
supports a non-zero energy current. There are Goldstone modes corresponding to this 
breaking of translational symmetry, and for D = 1 these actually destroy the long range 
order. The phase transition is then induced by the condensation of vortices, just like in 
the D = 1 planar rotor model. The Goldstone modes lead to a power-law modulation 
of the amplitude of the oscillatory spin correlation function. Even for D = 1, this state 
with only quasi- long range order is still capable of supporting a classical current. 

However, our model differs from those considered in in that, due to quantum 

effects, the current is also nonvanishing in the region < A < A c . This may be simply 
understood as follows. In @,0 the space integral J of the current happens to commute 
with the hamiltonian H^°\ This means, in particular, that the vacuum state of is 
also an eigenstate of H = — XJ. In the general case where the spectrum of 
has a gap, we expect the vacuum state of H to depend analytically on A, and therefore 
to remain unchanged at least for some finite neighbourhood of A = 0. In this region, 
the current will continue to vanish, and, indeed, all the equal-time correlations will be 
identical to the case when A = 0. The onset of the current then occurs at the first 
nonanalyticity, A = A c . In general, and in our model in particular, there is no a priori 
reason for and J to commute, and the structure of the vacuum may change as soon 
as A 7^ 0, allowing the current to flow. Interesting enough, we find that this happens 
due to quantum processes which first arise only at 0(h 3 ), so that such effects might be 
very small in a real system. 

Most of our analysis applies to the case when A is switched on in the disordered, 
symmetric phase of the continuous spin Ising model. However, a very similar picture 
seems to hold in the ferromagnetic phase. There is a transition at a finite value of A to 
a non-translational invariant state (for D > 1) which carries a classical energy current. 
However, in our particular model, in weak coupling, it can be shown that the expected 
second order transition is preempted by a first order one to a modulated state in which 
the higher harmonics are important. 

The outline of this paper is as follows. In the next section, after setting up an 
appropriate hamiltonian which captures the physics we want to describe, we calculate 
the effective action at tree level. Sec. 3. discusses the theory and its quantum corrections 
in the ferromagnetically disordered phase, first at small values of the bias, then at the 
critical value, and finally in the modulated phase. In all three cases we evaluate the 
classical current and its lowest order quantum corrections. In Sec. 4 we pay particular 
attention to the Goldstone modes which appear in the modulated phase for D > 1 and 
discuss how these dramatically modify the physics for D — 1. The next section describes 
our (limited) results for the ferromagnetically ordered phase, where there appears to be 
a first-order transition to the modulated phase at mean-field level. Finally, in Sec. 6 we 
summarise our results, contrast them with those found in the integrable D = 1 example 
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of Ref. 0, and connect them with those in other different, but mathematically similar, 
models of modulated phases. 

2 A lattice model 

Our first aim is to select a suitable field theory hamiltonian whose long-distance 
behaviour is expected to be in the same universality class as the transverse Ising model 
whose D — 1 version was discussed in ||. In the absence of any Lagrange multiplier, it 
is well known that the continuum relativistic 4 field theory, with a suitable momentum 
cut-off, has this property, and that it undergoes a zero-temperature phase transition in 
the same universality class as the classical Ising model in D + 1 space dimensions. 

However, it turns out that such a relativistic continuum theory cannot exhibit the 
physics we are trying to study at non-zero A. This is because the energy current is given 
by the components T Qi of the energy-momentum tensor, which, in a relativistic theory, 
are equal to those of the momentum density T iQ . The current J is therefore nothing but 
the total momentum P, which commutes with H^°K Any simultaneous eigenstate has, 
for large P, an energy E ~ P. From the convexity of the relativistic dispersion relation 
it is easy to see that the ground state of — \P is the usual vacuum for A < 1, while 
for A > 1 this operator is unbounded from below, so that the theory does not exist. 

Clearly it is necessary to modify the dispersion relation, and the natural way to do 
this is to consider the same field theory on a lattice, since this also provides an ultraviolet 
regulator. The hamiltonian of the D-dimensional lattice model that we consider is thus 

^ (O) = ^E{^ + ^ + E[0r-0r-eJ 2 + ^} , (1) 

where summation is taken over the vertices of a .D-dimensional cubic lattice and e a 
are unit (lattice) vectors. The lattice constant is taken to be unity. The canonical 
momentum, II, is introduced and the time coordinate is kept continuous to facilitate 
transformation to the lagrangian formalism later. 

Let us calculate the energy current jr a °f this hamiltonian. This satisfies 

z[^°),4 0) ] = V Q ^ Q , (2) 

where V a denotes the lattice divergence. The solution of (fj) for the space integral of 
the energy current, the only quantity of interest for us, is unique. 
After a simple calculation we find from 

Jr, a = \ [(n r + n r _ e J(0 r - fo-e*)] • (3) 

Note that this is hermitian but is odd under time-reversal and parity. 

Now we modify our lattice hamiltonian by the addition of the space integral of a 
component of the energy current, J2r 3r i> multiplied by a Lagrange multiplier, to get 

H=\T,{K + + £ [<h - 0r-ej 2 + ^4 + £#1} , (4) 
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Note that we have selected the first coordinate as the direction of the energy current. 

The term Y^rJri does not commute with H^°\ due to the presence of the 4 inter- 
action term. In fact the commutator is proportional to 

H^r^r+d -0r-ej (5) 
r 

Note that in the naive continuum limit this is proportional to 3 <9i0, which is a to- 
tal derivative, but when higher order derivatives are taken into account it has a non- 
vanishing effect. 

As we emphasised in the Introduction, the main advantage of our field theoretic model 
is that we are able to utilise the lagrangian formalism. However, it is instructive first to 
examine the hamiltonian @ in the limit g = when the self-interaction is neglected. In 
that case H may be diagonalised in terms of momentum-space annihilation and creation 
operators, with the result 

H = J^(^j:(l- C osp a )+^ 2 -Xsmp i y\p)a(p), (6) 

where the integral is over the first Brillouin zone. For fi 2 > and small A, the effect is 
only to tilt slightly the usual one-particle dispersion relation (solid line), as illustrated by 
the dotted line of Fig. ^. The ground state is however still the original vacuum. However, 
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Figure 2: Dispersion relations at A = 0, A = 1 < A c , A = 2.2 > A c (solid lines), and for 
the 'negative energy branch' at A = 2.2 (dashed line). 

there is a critical value of A > for which the dispersion curve touches the zero-energy 
axis at some positive pi = p^\ and beyond this value there are negative energy single 
particle states (dashed line of Fig. §), which will carry a nonzero energy current. Unlike 
the case considered in [[5] , however, where the elementary excitations were fermions, the 
particles will now try to Bose condense into the lowest energy single particle state. Since 
the number of such quanta is unbounded, the energy will diverge to minus infinity, an 
unphysical result. Of course, this picture will change in the presence of the repulsive 4 
interaction, since this will limit the total number of quanta. However, the interaction will 
also introduce scattering between quasiparticles of different momenta, and in particular 
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will introduce virtual processes involving the branch of the dispersion relation found by 
taking the other sign of the square root, indicated by the broken line in Fig. |[ For 
A greater than its critical value the energies of some of these states begin to overlap 
with those of the original single-particle states, indicating a complete breakdown of the 
quasi-particle picture. 

For this reason, among others, we have found it more useful to analyse the interacting 
theory through the path integral and the effective action formalism. First we need to 
calculate the lagrangian. We use the canonical equation to find 



Then the lagrangian, 



L = £n r 0r.-#, (8) 
r 



takes the form 



- Z)[0r - <Pr-e a f - JZ<t>rj 



(9) 



In the lowest order of the loop expansion the effective action is identical to the 
classical action. Then, at tree level, the classical field minimises the classical action. As 
we have argued above, for A larger than a critical value we expect to find a macroscopic 
number of particles with momentum pa in the ground state. This means that the 
classical field which minimises the action will be space-dependent. From the point of 
view of the path integral, there is also no reason to exclude a time-dependence as well. 
We shall assume that, since we are considering a stationary process, that this is not the 
case. Further justification for this may be found by considering the euclidean action, 
whose minimum is guaranteed to give the lowest energy state. This is found from (|9|) by 
changing the sign of all the terms and also letting — > icf). In general, then, euclidean 
saddle point solutions will be complex, unless they are time- independent. Since we are 
seeking a ground state with a real field, we should exclude complex solutions. 

If we allowed the classical field to depend on time the solution of the extremisation 
problem would become far from unique. In fact, one could find a solution at all values 
of the parameters, A and fi 2 . Furthermore, on a more practical level, we would com- 
pletely lose the guiding principle of finding the classical field by minimising the effective 
action. One can see this problem arising already at the classical level. If the classical 
field would depend on time, then its Fourier transform would depend on 'energy' and 
the kinetic part of the lagrangian would become indefinite. Looking at this question 
from a slightly different angle, remember that the second functional derivative of the 
effective action, or at tree level, classical action, is the inverse propagator. The inverse 
of this, the propagator has a spectral representation with a positive definite weight func- 
tion and non-negative threshold. Thus, at any space or timelike four-momentum (or 
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generalised four momentum for lattices) the propagator is positive. In momentum rep- 
resentation this insures the definiteness of the inverse propagator, as well. However, if 
the energy component is nonvanishing then the spectral representation does not insure 
the definiteness of the propagator and there is no minimum principle for the effective 
action. 

Note that, as we shall discuss in detail later, the fact that the classical field is time- 
independent does not preclude the ground state from carrying a current. This is because, 
as may be seen from (0,0), if <fi is space-dependent then <fr may vanish without IT, and 
therefore the classical current, vanishing. From the hamiltonian perspective, this is re- 
lated to the fact that the 'time' appearing in the lagrangian (§) is not in fact the true 
time, but rather the conjugate quantity to the effective hamiltonian H = — \J, 
whose ground state properties we are trying to find. The real time-dependence of the sys- 
tem is still generated by the original hamiltonian H^ '. Thus, in the Heisenberg picture, 
having a 'time'-independent classical field which commutes with H in fact corresponds 
to a true time dependence = — i\[J, 0] - thus, as expected, local inhomogeneities in 
are transported along by the current in this state. 



3 Energy current in the disordered phase fj 2 > 0. 

The minimisation of the classical action?] on functions independent of time leads to the 
following equation for classical field 0^1 that will be henceforth denoted by <pr'- 

.5T M A 2 



i ~Fr' = ~T l 2( ^ r ~ ( Pr-2e 1 ~ <Pr+2e 1 \ 
o<pr 4 

D 



- E[20r--0r--e Q -0r+eJ-^0r-f[0r] 3 = O. (10) 

Introduce the Fourier decomposition of the time dependent classical field 0r 

<Pr = (2?r)(D+1)/2 J ^ d D p J Q dE ty(p, EyP- r ~ lEt + h. c] . (11) 

and let <p{p, E) = ip(p, E) + ^{-p, -£).[] 

A time-dependent solution for the classical field then has the form setting <f)(p, E) = 
27r(p(p)S(E) ( |T0"[ ) where 0(p) satisfies 

= ° = D o _1 (P,0)^(P) " f [<^0p)> (12) 
and we have introduced the notation 

D 

D \p, E) = (E + \ sm Pl ) 2 - 2 ]T(1 - cos Pa ) - /i 2 (13) 



3 For convenience, we will use the negative of the effective action that will have a maximum, rather 
that a minimum at the correct classical field. 

4 This implies the relations 4>{p, E) = ip(p, E) for E > and <f)(p, E) = tp^(— p, —E) if E < 0. 
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and 



k ](p) = J d D p 1 ...d D p k 5(p-p 1 - (14) 



(27T) 

The tree level expression for the inverse propagator is 
Z 6<f>(p,E)6(f>(-p>,-E>) 



5(E - E') [s(p - p')D,\p, E) - |[0 2 ](p - p')} • (15) 



Note that i?) vanishes on the single particle dispersion curve discussed earlier, 

as well as on the 'negative-energy' branch. The coefficient Dq X (p, 0) in (|l"2"D can be easily 
analysed. It may have several extrema. Besides the one at p = 0, provided that A > 1, 
there are two at p\ = ±Pi°\ p a = 0, if a ^ 1, where cospf^ = 1/A 2 . At A > 1 the latter 
ones are the only maxima. The value of Dq 1 at these maxima is 

£oV 0) ) = (A-l/A) 2 -^ 2 . (16) 

Then for every real value of the mass gap, fi, there is a A c > 1, such that the maximum 
of this coefficient is exactly zero, 

1 



Note that the convexity of the effective action requires that after setting E = 
(|j~5|) is negative at the appropriate classical field. Indeed, when Dq X (p) > one can 



find a nonvanishing solution of ([12]). In fact, this nonvanishing solution will provide 
the minimum of the action in contrast to the symmetric solution, = 0. We expect 
that if the vacuum expectation value is independent of time then the propagator (the 
second derivative of the action) is negative definite. Then the effective action is a convex 
function of the classical field and the classical field should be chosen to minimise it. 



3.1 The symmetric phase 

Suppose first that A < A c . Then Dq 1 ^ ^) < and the only real solution of ( |I2| ) is 
<f)p = 0. Consequently, the symmetry is unbroken. The tree level inverse propagator 
becomes 

D \p,E) = Ao\p,E + \sm Pl ), (18) 

where A is the propagator in the absence of the perturbation. Then the only change 
in the perturbed theory is the change of the energy-momentum dispersion relation to 



E = — A sinpx + m, 2 + 2 ^(1 — cosp Q ). 

Though this relation dips near p = p(°\ the energy is always positive. In other words, 
the energy gap is nonvanishing. Perturbative corrections are not expected to change this 
result, but they may change the location of the critical point. 



S 



Although the propagator is modified when A 7^ 0, in fact the effect on the ground 
state properties is rather subtle. For example, the equal-time spin correlation function 
in the absence of loop corrections entails integrating Dq(p, E) over all E, and therefore 
the correction term Asinpi is simply shifted away. Does this result survive the addition 
of loop corrections? This may be examined by making a time-dependent transformation 
on the Fourier components of the field 

<P(p,t)^<p(p,t)e UXsi ^\ (19) 

which shifts away the Asinpi term in the bare propagator. The interaction now has the 
form 

4 4 
J rftn[4 (n V(P (n) ,t)]e lUE »= lSm!, '" 1 (5(Ep (n) ). (20) 

n=l rt=l 

From this it can be seen that if we approximate sinp! by pi, the t-dependence disappears 
and the result is the same as in the theory with A = 0. This is clearly connected with the 
fact that, to lowest order in the derivative expansion, the current does commute with 
the interaction. 

The subtlety of this modification may be seen, for example, in the computation of 
the expectation value of the current in this phase. From (^, [7]) this is 

(if) = J dEd D p{E + Xsm Pl )X sin p^p, E^-p, -E)) . (21) 

If we make the above shift E — > E + Asinpi, and we compute the correlator in the 
absence of any interaction, it is easy to see that this vanishes through the symmetry 
E — > —E (or pi — ► —pi). Indeed, this will remain true even when the interaction is 



included, if we can ignore the momentum-dependent phase factors in (|20[ ), since these 
symmetries then continue to hold. But the phase factors couple these two symmetries 
together, and, since the current is even under the joint reflection of E and p\ it does 
not necessarily vanish when loop corrections are included. However, it turns out that 
the first non-zero contribution first occurs at three loops, due to the diagram shown in 
Fig. |3]a. The one- and two-loop diagrams still vanish. For example, at the internal vertex 
in the diagram shown in Fig. |3|b the incoming momenta are (p, — p,p', — p') so the phase 
X^=isinp n still vanishes. We conclude that the first non- vanishing contribution to the 
current in this phase occurs at three loops, and so is 0(g 2 h 3 ). 

3.2 The critical theory 

At the critical point A = A c , the classical field still vanishes at the minimum of the 
action. Therefore the tree level propagator has the form 

D<yP ' E) = (E + X c sin Pl ) 2 - 2 £ a (l - cosj9 a ) - p 2 ^ 

It is possible to make the same transformation as above and shift E in the bare propa- 
gator. However, this obscures the fact that there is now no energy gap, and the theory 
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a 



b 



Figure 3: Diagrams contributing to the expectation value of the energy current: a) Lowest 
order nonvanishing diagram, b) Vanishing two loop diagram. The cross corresponds to 
a bare energy current insertion. 



may develop infrared singularities close to E = which modify the large distance be- 
haviour of the correlation functions. This can be seen if we expand the denominator of 
the propagator around p = ±p^ and E = 0. The constant term vanishes at the critical 
point. The term linear in p also vanishes, as ±p(°> were defined as locations of maxima 
of the inverse propagator. Then, keeping leading order terms in E and p =F p^ only, we 
obtain 

D(P, E) ~ — ; m ; 2 l m , J^-^ (23) 



±2£A c sin P r - \1&?pV(Pi TpD 2 ~ E^iti 



pf = cos- 1 1 = cos- 1 ( 1 - p^-Hz£ ) > 0. (24) 



where 



The two signs in fl23|) correspond to poles that are separated from in other by a finite 
distance. These two poles lead to two propagators of states which, at least in the absence 
of interactions, do not mix. After rescaling E, as 2EX C \ sinp[ | — > E, and shifting and 
rescaling pi as A c sinpJ (pi — Pi) — > pi we obtain the following expression, which can 
be regarded as propagators for a 'particle' A and for an 'antiparticle' A, both, however, 
non-relativistic: 

When the interactions are included, in terms of the non-relativistic excitations above, 
the terms allowed by overall momentum conservation are the scattering processes AA — > 
AA, AA — ► AA, and also AAAA — > and its inverse. Without these last two processes, 
the theory is identical to that of non-relativistic particles with a local 2-body repulsive 
interaction. Such as theory does have infrared singularities below D = 2, but they do 
not affect the propagator, since it is easy to see that there are in fact no loop corrections. 
For the same reason the quantum corrections to the current all vanish. 

This is no longer true in the presence of the last two vacuum processes, which indeed 
lead, at two loops, to corrections to the current corresponding to those illustrated, for 
fi 2 > 0, in Fig. [3]a. However, because of the way the momentum flows through these 
diagrams, it can be seen that they are not in fact infrared singular. Thus we expect 
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that the quantum corrections to the current are not singular, leaving the singularity in 
classical current, to be studied in the next section. 



3.3 The broken phase 

When A is increased beyond its critical value, A c , then the maximum of the inverse 
propagator becomes positive, D ( ^ 1 (p < - - ) ) > 0. Then a symmetry breaking solution of 
(|l2l) exists, along with the symmetric solution, <f)(p) = 0. As the nonzero solution 
minimises the effective action we must choose it over the symmetric one. 

It is not a priori clear whether the relevant solution of (|l|) is a classical field with a 
sharp spectrum (i.e. delta function in momentum space) or at any given A the spectrum 
of the classical field contains all the momenta for which D~ l (p) > 0. Examining the 
second functional derivative resolves this problem. On one hand, the convexity of the 
action demands that the right hand side of (|T5|) is negative even in the neighbourhood 
of p = p^ . On the other hand, at A > A c the multiplier of the momentum conservation 
delta function (that is D~ 1 (p)) is positive in the neighbourhood of p^°\ The negativity 
of the second variation can only be secured if the last term of ([12]) also contains a 
term proportional to 5(p — p'). In turn, that is possible only if <j)(p) itself contains a 
term proportional to a delta function (sharp spectrum), as well. Suppose that the delta 
function fixes the momentum at p = ko, or in other words, ~ 5(p — ko). Then it is 
easy to see that the minimum of the action, S m i n ~ — [D~ l {ko)} 2 - This is minimised, as 
a function of ko if ko = ±p(°\ In other words, the delta function in the classical field 
must set the momentum equal to ±p(°). 

Since <f)(x, t) is real at E — both "negative" and "positive" frequency terms con- 
tribute. Therefore, we use a symmetric ansatz 



In (|27|) p is a real amplitude, and 9o is an arbitrary real phase. As we shall see in 
the next section the arbitrary phase is related to the presence of a Goldstone boson in 
the broken phase. 

Near the phase transition the term A0(p) is much smaller than the leading term, 
4>o(p)- The self-consistency of this requirement will become obvious below. 
Upon substituting fl2"6"|) into condition ( P"2"|) we obtain 



4>{p) = <Mp) + A0( p ), 



(26) 



where 



o (p) = p \ e te ° 5(p - p (0) ) + e- ie ° 8(p + p {0) ) 



(27) 








+ A 2 sin 2 p x - 2 



cosp a ) - p 2 A0(p) 



a 




3p (o)) + e - 3ie °5(p + 3p (0) )l + 0(p 2 A(f)(p)). (28) 
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Consequently, the leading order term (0(X — A c )), proportional to 5(p±p^°>) of (|2 
vanishes if we choose 

X x2 



p2 = 2(2vr; 



D 

(29) 



A--) 
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Then near the phase transition the amplitude of the classical field vanishes as p 
\/A — A c . Note that at p 2 = A c = 1 and p ~ A — A c . 

The leading order third harmonics terms also cancel if we choose 



A0(p) = p 3 [e My °6(p - 3p w ) + e-^°5(p + 3p w )J + 0(p 5 ), (30) 
where the amplitude of the third harmonics, p^, has the form 

Ps = ~ 2(1 - cos 3p(°)) - A 2 sin 2 3p(°) + p 2 Q{2tc) d P ' ^ 

It follows from (§T]) that A0(p) = 0((A - A ) 3/2 ). 

In general, it is easy to see that near the phase transition A ~ A c one can solve 
equation (|T2"D iteratively, resulting in amplitudes for the higher harmonics that have a 
following behaviour near A = A c : 

fe ,«~(A-A(»>) (2 " +1,/2 . (32) 

Here p2 n +i is the coefficient of (2n + l)st harmonics of momentum p(°\ 

The conclusion is that, at least at the classical level, the system undergoes a second 

order phase transition to a phase with spontaneously broken translation symmetry. 
Finally, we comment on the classical part of the energy current, found by substituting 

the classical expectation value of the field into (^). Substituting the canonical momentum 

by its lagrangian expression and using that the classical field is time independent we 

obtain 

( J d D xj*(x)) = -y / d D x^ x+ei ) - fo-cJ (te+c) - <Px-e a ))- (33) 

This expression is obviously zero unless a — 1. For a = 1 we obtain, after substituting 
(f) x = p ^- x + h.c. 

</^»> * VAWp<°> ~ v ~ \ { (A A : A A > l$H (34) 

As well as this 0(1/ g) contribution, there are perturbative quantum corrections just as 
for A < X r . 



4 Goldstone modes and the case D = 1. 

Up to this point we have investigated the D > 1 case only In one dimension, the 
Mermin-Wagner-Hohenberg theorem asserts that there can be no spontaneous breaking 
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of the continuous symmetry of translational invariance, due to infrared singularities of 
the Goldstone modes. In fact, the physics is very similar to that of spin or charge 
density waves, and of modulated phases such as occur in the ANNNI model 0, in 1 + 1 
dimensions. 

Rather than include all the fluctuations about the classical solution, it should be 
adequate to derive an effective action for these Goldstone modes only, by writing the 
field <j)(x, t) in the following form: 

<t>{x,t)=p{[l + t{x,t)}e i9 ^e ipi ° )x }+c.c., (35) 

where p and are chosen to minimise the classical action, and £(x,t) and 6(x,t) are 
real. These functions are chosen so that their fourier transforms have support only for 
\p\ < p(°>: in this way the decomposition in fl35| ) is unique. In principle we should 
also include higher harmonics, but, as shown in Sec. 3.3, these are suppressed near the 
transition. We have included possible longitudinal fluctuations through the field £(x,t). 

Upon substituting (^) into the action and keeping terms quadratic in oscillating 
fields £(x,t) and 9(x,t) only we obtain 

s = whlJ dDpdE 

{[D \p - p®,E) + D \p + p(°\E) - 3gp 2 } £(p, E) £*(p, E) 
+ [D \P ~ P {0 \ E) + D \p + p(°), E) - gp 2 } 6(p, E) 6*(p, E) 
+ i [D \p - p(°),E) - D \p + p(°\E)] [6(p, E) C( P , E) - 6*(p, E) £(p, £?)]}(36) 

where the notation Dq X (p, E) — (E + X smpi) 2 — 2 ~ cosp a ) — p 2 has been intro- 
duced. Terms containing fields £(p,E) or 9(p,E) taken at a momentum of 0(p^) were 
omitted from (|3~6|) . For convenience, we also symmetrised (|36|) in momentum. 
Expression ( j36|) of the action has a few important features: 

• p measures the deviation of the momentum from its "critical" values, ±p^ . 

• Radial and angular degrees of freedom are mixed. The term mixing them is pro- 
portional to A. Thus, no such term appears in the somewhat similar case of spon- 
taneous breaking of U(l) global symmetry, where the radial and phase degrees do 
not mix at tree level. 

• If we omit the mixing term then, as expected, the radial mode is massive, while 
the angular mode is massless, In other words, the energy gap vanishes at p = 0: 
D - 1 (p(°),0) + J D 1 (-p( ),0)-^ 2 = 0. 

Thus, due to the mixing of the two modes, which is a peculiar feature of our model, we 
cannot analyse the zero modes by omitting all degrees of freedom but the phase. The 
eigenmodes of propagation are found by diagonalising the propagator matrix in the two 
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dimensional space (£,#). The eigenvalues of this matrix are 



2 

£>-i = L ^D \p - p {0 \E) + D^\p + p (0) , E) -2gp 2 ± \/g 2 p 4 + 16X 2 E 2 sinV 0) cos 2 p} 

(37) 

In the far infrared region, E,p 2 << gp 2 , one finds the mode corresponding to the upper 
sign tends to the #-mode, (the £ admixture is of 0(E/gp 2 )), while the one corresponding 
to the lower sign tends to the £ mode. Yet the dispersion relation obtained for the upper 
sign is different from the one we would have obtained, had we omitted the "radial" mode 
from the outset. This can be seen if we consider the asymptotic form of eigenvalues of 
the inverse propagator matrix, 

D . l= ( (v-im(l&- P w)-p'pl \ 

1,2 \ [p 2 - I (A 2 - 1/A 2 )] E 2 - (A 2 - 1/AVp? - p 2 p 2 - 2gp 2 ) 1 ' 

The first of these modes is massless. The second mode appears to be massive, except 
near the phase transition, if p 2 < 4(A 2 — l/X 2 )/g, the energy squared term has the 
wrong sign. Having imaginary energy, this is not a propagating mode. We will soon 
show, however, that at small values of p the classical approach becomes unreliable, at 
least at D — 1. 

Next we shall consider the case D — 1. Near the phase transition the regularised 
equal time correlation function of the gapless mode becomes 

D( - \ 1 f rlP rl eiP(X ~ V) ~ 1 y/9log\x-y\ /qqA 

[X V) ~ i47r 2 (A 2 -l/A 2 )i P AE 2 /g-p 2 p 2 ~ 8vrp(A 2 - 1/A 2 ) 1 ' 

The correlation function is proportional to 1/p, so it diverges as 1/ y/X — X c near A = A c . 

Near E — (|5D|) approximates the correlation function for the phase mode of oscil- 
lations. Then the equal time correlation function of <p(x, 0) is 

Mx, O)0(i/, 0)) ~ ^M-^M+^M) ~ p 2 cos[p^(x -y)} ^ 

\x — yyi 

where r] = ^fg / '[8tt p(X 2 — 1/A 2 )], which exhibits an oscillating power-like behaviour with 
a nonuniversal critical exponent. This means that there is only quasi- long range order 
in D = 1. 

Note that the parameter p oc (A — Ao) 1 ^ 2 plays a very similar role to (3 = 1/T 
in the XY-model. Just as in this model, there are also vortex-like excitations whose 
condensation at high enough T (small p) completely destroys even the quasi-long range 
order. In our case, these result from the fact that 8(x,t) is an angular variable, defined 
only modulo 2tt. A 'vortex' corresponds to a fluctuation in which an extra oscillation of 
the field <f)(x,t), localised in x, is either inserted or removed in a large but finite space 
and time interval. This has an action which grows like the logarithm of the system size. 
By the Kosterlitz-Thouless criterion ||10|| , they become relevant when rj > 1/4. 

We conclude that the transition to the quasi-modulated phase occurs not at A ~ A , 
but at a higher value when p has become sufficiently small. At the transition, we should 
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then have r\ = 1/4, with the value of this exponent decreasing as we go further into the 
modulated phase. 

Although this phase exhibits only quasi-long range order, it nevertheless supports a 
non-zero classical current, as for D > 1, since the leading contribution ~ Xp 2 sin 2 is 
not affected by infrared problems. In fact, this quantity corresponds to the superfluid 
density in the XY model, and is expected to jump discontinuously at the transition |JTT[ . 

In the ultraviolet limit of ([J7|) the two roots of the inverse propagator reduce to the 
two poles of fl25|). This shows the consistency of our descriptions of the critical and 
broken phases. 



5 Energy current in the ferromagnetically ordered 
state. 

In the regime p? < we have not been able to obtain results comparable to those above, 
especially in the difficult case of p 2 = 0. We can show however, at least at the classical 
level, that the second order transition at p 2 > turns into a first order transition at 
p 2 <0. 

When p 2 < the Z 2 symmetry of the theory is broken even at A <C 1, where the 
translation invariance is certainly not broken. Thus, we must shift the field to minimise 
the effective potential, but the classical field is a constant, (f) Q . When A is increased the 
classical field will become coordinate dependent above a certain critical value. Writing 
the classical field as 4>(Pi) = 4 ) o + i'i'Pi) the extremum condition for field ip will read as: 



5T 



eff 



Slf)(p) 







D 



A 2 sin 2 ]?! — 2VYl — cosp Q ) — 2\p 2 \ 

a 

2(2^p J d D pid D p 2 ip{p 1 )i){p 2 )5{p-p l -p 2 ) 

q \ r 

6(2^p J dDpidDp2dDp3 ^( p ^( p ^( p ^ 6 ( p ~ P± ~ p z -PsM 41 ) 



where 0o = \/6|A* 2 |(27r) I) /(7. 

We can rescale (|4lD by dividing by \p 2 \4>o and introducing the new field ip —> ^^Q- 
Then, using the fact that the ground state depends only on pi (in what follows we will 
drop the subscript) , one obtains the following equation for 







A 2 sin 2 p — 2(1 — cos p) 



1p(p) 



+ 3 / dp 1 dp 2 ip(p L )ip(p 2 )5(p — p 1 — p 2 



+ 



dp 1 dp 2 dp 3 ^p(p 1 )ip(p 2 )^/j(p 3 )S(p — p 1 — p 2 —p 3 )- 



(42) 



Naively, one would think that just like for p 2 > 0, at least at the tree level, the 
phase transition is determined by the linear term in ((42"|) and occurs at A = Ai, where 
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Ao — 1/Aq = \/2/x. (Note the extra factor of \/2, reflecting the fact that the bare mass for 
/i 2 < is a/2|/x|.) Contrary to this expectation we will prove that, due to the presence 



of the term quadratic in ip in (42), the system undergoes a first order phase transition, 
at a critical value of A = A c < Ai, to a state with broken translation invariance. In order 
to show this we merely have to exhibit a solution of (|42|) which gives a lower action, for 
A c < A < Ai, than the solution ip{p) = 0. To this end, consider the ansatz 



ipo8(pi) + pi 



+ ^2 



A2B 



e^Sfa. -Pi) +e 
- 2pJ 0) ) + e^fa 



-i0, 



2p\ 



(0)^ 



(43) 



where, as before, cosp^ = 1/A 2 . After substituting this ansatz into 
the action with respect to x/jq, pi, and ip2- We obtain 

-2^o = 3(^ 2 + 2p 2 + 2^ 2 2 )+^ 3 + 3^o(Pi+^i 2 ) 

epi = 6(^oPi +4>2pi) +3pi(pi + 2ipl) 
-2^ 2 = 6(p 2 + ^o) + 2^oP 2 + ^oV>2 



we minimise 



(44) 



where we used the relation 



cos(2p^°' ) )] = and the notation 



A 2 sin 2 (2pS 0) ) - 2[1 

p 2 | - 2 = e. 

Note that in terms of the new variables a second order transition is possible only at 
e > 0. (|44]) can however be solved numerically and doing so we find that it has a real 
solution for all e > e = —0.3395. This value of e corresponds to A c ~ 1.834 < Ai = 
1.932. Below A c = 1.834 there is no real solution of the form ([43"D , but above this value 
a pair of complex roots become real and the parameters take the values pi ~ 0.130, 
ipQ ~ —0.0410, and ip2 — —0.0358. The true tree level transition is then at a value 
A c < 1.834. The phase transition line is at a fixed value of (A — l/A) 2 /|p 2 |. This fixed 
value is smaller than 1.6605, and certainly smaller than the superficial prediction, 2, 
that was based on the assumption that the transition is second order. 

Of course, this argument does not imply that the ansatz (fl3"D represents the true 
minimum of the effective potential (now there is no reason to ignore higher harmonics 
as there was in the second order case) but is does imply that the transition cannot be 
second order, at least at tree level. Nor does the argument rule out the possibility of a 
restoration of the second order nature of the transition by the fluctuations, especially in 
low dimensions (see next section). 



6 Summary 

Following investigations of similar nature for integrable one dimensional spin chains || 0, 
we have studied the effect of an energy current in field theories, at zero temperature. 
The energy current was imposed on the system by adding the global energy current op- 
erator, multiplied by a Lagrange multiplier, A, to the Hamiltonian. The common feature 
of these models is the appearance of a phase transition at a critical value of A. 
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Field theoretic models modified by an energy current can roughly be divided into two 
groups. In the first, less interesting group of theories the hamiltonian is destabilised at 
the critical value of A, A c , whereby the energy spectrum has no lower bound at A > A c . 
A typical example for these models is a relativistic, self-interacting, real, scalar field 
theory. Theories, in which at D > 1 (where D is the number of spatial dimensions) 
translation invariance breaks spontaneously, are more interesting. A typical example for 
this second group is a self-interacting real scalar field theory on a lattice. 

In the lattice models that we consider the global energy current does not commute 
with the Hamiltonian and consequently, unlike in spin models [[J 0, the ground state 
expectation value of the energy current has a nonvanishing perturbative contribution at 
all A 7^ 0. Another, nonperturbative contribution to the energy current appears at the 
phase transition point. 

The system has a few unusual features in the broken phase, at D > 1. The first of 
these is that the field has a nonvanishing, oscillating vacuum expectation value. Cor- 
relation functions are also oscillating functions of the coordinate. The oscillations are 
the consequence of the unusual form of the energy-momentum dispersion relation and 
the vanishing of the energy gap at a nonvanishing value of the momentum. The second 
order phase transition to the broken phase is accompanied by the appearance of a Gold- 
stone boson, in spite of the fact that the translation symmetry is a discrete symmetry 
on the lattice. This happens because near the phase transition the lattice structure of 
the system becomes irrelevant. 

The second interesting feature of the field theoretic model is the apparent sponta- 
neous generation of degrees of freedom and enlargement of the internal symmetry near 
the phase transition. As the energy gap vanishes at p = p° ^ 0, it is meaningful to 
rewrite the real scalar field into the form <f>(x) = x(x)e ip( ° : x + x 1 (x)e~ tpm . At least in the 
infrared domain (p << p^) x{ x ) is a genuine complex field with U(l) symmetry (planar 
model) in contrast to the Z2 symmetry of the of the field <f>(x). When the vacuum ex- 
pectation value of x is non-vanishing, its phase degree of freedom becomes a Goldstone 
mode. The extra near infrared degrees of freedom can of course be related to large 
momentum modes of the original theory, so in reality the number of degrees of freedom 
is unchanged. What we really observe is the transmutation of a space-time symmetry 
(translation invariance) into a broken internal symmetry (phase transformations of field 

x). 

The analogy with the U(l) model in D > 1 spatial dimensions is not complete. Due 
to the peculiar form of the single particle dispersion relation the propagator takes a 
nonrelativistic form in the infrared domain. This leads to a critical dynamics different 
than that of a complex scalar field, and the anomalous dimension vanishes at all D > 1. 

When D — 1 the translation invariance cannot be broken spontaneously. Yet the 
characteristic momentum p^> and the complexification of the degrees of freedom survive 
due to large fluctuations, the field x( x ) m ust have a vanishing position-dependent vac- 
uum expectation value. Still, the analogy with the U(l) symmetric XY model, based 
on the decomposition of the scalar field into a superposition of a complex field and of its 
complex conjugate, is maintained. One obtains damped oscillating correlation functions 
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with non-universal anomalous dimensions. The phase transition at A = A c is driven by 
the condensation of vortices of the x( x ) field. 

Despite some similarities with the exactly solvable cases treated in || @, there are 
also some important differences in our results. Partly this is due to the fact that in our 
case the current does not commute with the hamiltonian, so that the current is non- 
vanishing even for arbitrarily small bias A. However, the singular part of the current, 
which we have argued occurs only in its 'classical' piece, has a different behaviour, 
vanishing linearly as A — > A c +, rather than proportional to \/X — X c as in 0. In D = 1, 
this linear vanishing is modified by the vortices to a discontinuous jump. This behaviour 
is related to that of the correlation function: we find that in our model it should exhibit 
oscillations modulated by a power law | j 77 , where r\ is continuously varying throughout 
the modulated phase, increasing to the universal value 1/4 at the transition. This should 
be compared with the constant value rj = 1/2 found for the transverse Ising chain in ||. 

These differences in detail are not necessarily surprising. When such integrable mod- 
els as the Ising or XY chains are mapped onto coulomb gas models, the additional 
conservation laws in the integrable system may show up in unusual constraints (for ex- 
ample the charge at infinity) which can modify the critical exponents away from their 
usual values. However it should be clear that our results do not depend so much on 
details of our model, and therefore should possibly be more generic. 

A comment should be made concerning the time-dependence in our model. In most 
of this paper the parameter t referred to as 'time' is merely an artifice introduced in order 
to provide a lagrangian formulation of the problem of finding the ground state properties 
of the operator — A J. Thus no meaning should be attached to quantities which are 
not t-independent. The true time dependence is still generated by H^°\ In particular, 
this means that the expectation value of the current, which does not commute with H^°\ 
should be time-dependent. The constant quantity which we have computed is therefore 
only its zero-frequency part. It would be interesting, but difficult, to investigate the 
time-dependent part as well. 

In fact, the mathematics of our model bears a very close resemblance to that which 
appears in the study of other systems which exhibit spatially modulated phases, such as 
spin or charge density waves, or systems such as the ANNNI model [[J with competing 
ferromagnetic nearest and antiferromagnetic next-nearest neighbour interactions. In 
fact, if one neglects the 'time'-dependence in (|j), and takes the limit g — > oo, yU 2 — > 



— oo with fi 2 /g fixed, so that the field at each site is frozen to the values ±y3|/i 2 |/g, 
this becomes precisely the energy function for the classical ANNNI model in D + 1 
dimensions. Even in mean field theory the phase diagram for this model is extremely 

rich. Between the ferromagnetic •••++++••• phase and the • • •+- 1 !-+••• antiphase 

lie an infinite number of other modulated phases. In 1+1 dimensions, these can also 
be separated by Kosterlitz-Thouless incommensurate phases, bounded by Pokrovsky- 
Talapov transitions [12]. It should be stressed that the limit above is quite different 
from the regime considered in the body of this paper where it was assumed that the 
characteristic wave vector is much smaller than the reciprocal lattice vector, and 
that lattice effects are irrelevant. The modulated phases considered there are therefore 
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incommensurate. The 'time '-dependent fluctuations in our model are rather different 
from those in transverse direction of the ANNNI-type models. Nevertheless it is possible 
that for D = 1 the first-order transition predicted in the last section is replaced by one 
of Pokrovsky-Talapov type. 

Besides these correspondences, there remains the important question of the relevance 
of these kinds of model, whether integrable as in Refs. or non-integrable as in 

our example, to the study of non-equilibrium processes at finite temperature where 
dissipation plays a crucial role. It is to be hoped that some of the phenomena discussed 
in this paper will survive into this regime. 

Acknowledgements 

The authors wish to thank T. Antal, F. Essler, Z. Racz and J. Yeomans for discussions. 
J.C was supported in part by the Engineering and Physical Sciences Research Council 
under Grant GR/J78327 while P.S. by the U.S. Department of Energy under grant 
number DE-FG02-84ER-40153. 

References 

[1] M.C. Cross and P.C. Hohenberg, Rev. Mod. Phys. 65 (1994) 851. 

[2] B. Schmittmann and R. P. K. Zia, in Phase Transitions and Critical Phenomena, Eds. C. 
Domb and J.L. Lebowitz (Academic Press, New York, 1995). 

[3] S. Katz, J.L. Lebowitz, and H. Spohn, J. Stat. Phys. 34 (1984) 497; A. DeMasi, PA. 
Ferrari, and J.L. Lebowitz, Phys. Rev. Lett. 55 (1985) 1947. 

[4] T.L. Hill, J. Chem. Phys. 76 (1982) 1122; P.L. Garrido, A. Labarta, and J. Marro, J. 
Stat. Phys. 49 (1987) 551. 

[5] Antal, Z. Racz, and Sasvari, Phys. Rev. Lett. 78 (1997) 167. 

[6] T. Antal, Z. Racz, A. Rakos, and G. M. Schiitz, Phys.Rev. E 57 (1998) 5184. 

[7] G. Jona-Lasinio, Nuovo Cim. 34A (1964) 1790. 

[8] M. Peskin and Schroeder, "An Introduction to Quantum Field Theory," Addison- Wesley 
(Reading, Mass. 1995). 

[9] See, for example, J. Yeomans, Solid State Physics 41 (1988) 151; W. Selke, Phys. Rep. 
170 (1988) 213. 

[10] J. M. Kosterlitz and D. J. Thouless, J. Phys. C 5 (1972) L124; J. Phys. C 6 (1973) 1181. 
[11] D. R. Nelson and J. M. Kosterlitz, Phys. Rev. Lett. 39 (1977) 1201. 
[12] V. L. Pokrovsky and A. L. Talapov, Phys. Rev. Lett. 42 (1979) 66. 



19 



